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1 Introduction 



It has recently been claimed by Damour, Deser and McCarthy |Q] that the Nonsymmetric Gravita- 
tional Theory (NGT) || is unphysical due to radiative instability. They claim that NGT predicts 
an infinite flux of negative energy from gravitational waves. We prove this claim to be false by 
demonstrating that the flux of gravitational radiation from an isolated source in NGT is identical 
to that found in General Relativity (GR). 

Our proof is based on the following analysis. Firstly, we solve the NGT field equations for an 
isolated axi-symmetric source using the expansion technique employed in |TJ. (The exact version 
of this work will be published elsewhere |J). We find that the skew sector of the theory does not 
contribute to the flux of gravitational radiation due to the short-range nature of the skew fields. 
Secondly, since the NGT wave equations allow linear superposition in the wave zone, we see that any 
compact source has a wave field that can be expressed as the superposition of a suitable collection 
of differently aligned axi-symmetric solutions. This completes our proof. 

2 NGT Vacuum Field Equations 

The NGT Lagrangian without sources takes the form: 

c = V=gg^R»u(w), (1) 

with g the determinant of g^ v . The NGT Ricci tensor is defined as: 

R,u(W) = - + W^J - W&Wfr + W%W^ (2) 

where W^ u is an unconstrained nonsymmetric connection : 

(Throughout this paper parentheses and square brackets enclosing indices stand for symmetrization 
and antisymmetrization, respectively.) The contravariant nonsymmetric tensor g^ v is defined in 
terms of the equation: 

g^g*v = tf»g va = e (4) 

If we define the torsion vector as: 
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then the connection r^,, where 



is torsion free: 



rj„ = + -s*w v , (6) 



r„ = rf H = o. (7) 



Defining now: 



we can write: 



Rfiu(T) — (^'(u0),u + r fi//3),u) r a^ r u/3 + r f a m r uiy) ( 8 ) 



flA«/(w) = fl^(r) + |w M , (9) 

where Ifuj = \{W^^ — W v ^). Finally, the NGT vacuum field equations can be expressed as: 

9[±v,cr — gpb^^a ~ 9pp^au = (10a) 
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[H )„ = 0, (10b) 



R^(T) = \w {vM . (10c) 
It is convenient to decompose into standard symmetric and antisymmetric parts: R(w,) > R\u,v] > 



and then rewrite the field equation (10c) in the following form: 

iV)(r) = o, (Ha) 

%H,P}( r ) = %H,p( r ) + ^M,p( r ) + r \mA v ) = • ( llb ) 
We shall base our analysis on the expansion used in , where the field equations are expanded 
in powers of the antisymmetric part, h^j, of the metric but where the symmetric part is taken to 
be an exact GR background. The field equations become to lowest order: 

GR Rpu = 0, (12a) 
D a F tlua + 4 GR R a ^ u h [ ^ ] = ^W M , (12b) 
D»h [H = 0, (12c) 

where F^ ua is the cyclic curl of fy^] and is the GR background covariant derivative. 

For our current purposes, we shall be using the radiative, axi-symmetric GR background found 
by Bondi, van der Burg and Metzner 0]. Since their result is given as an expansion in inverse powers 
of the radial coordinate, we shall also have to expand the above equations in the same way. The two 
expansions are perfectly compatible, and both expansions are well suited to studying the equations 
in the wave zone. 
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3 The Metric 



(13) 



The physical situation that we shall be studying is that of an isolated spherical body which has 
been deformed by a perturbing compression along a single axis. We seek to study the flux of 
gravitational radiation emitted as the body oscillates through the cycle oblate-spherical-prolate- 
spherical-oblate. 

Due to the physical picture sketched above, and to the fact that we are interested in the 
asymptotic behaviour of the field at spatial infinity (in an arbitrary direction from our isolated 
source), polar coordinates x° = u, x = (r, 9, <f>) are the natural choice. The "retarded time" u = r — t 
has the property that the hypersurfaces with u = constant are light-like. Detailed discussion of 
the coordinate systems permissible for investigation of outgoing gravitational waves from isolated 
systems can be found in j|, ||]. 

The covariant GR metric tensor corresponding to the situation described above: 

/ Vr -i e W _ u 2 r 2 e 2 ^ e 2 ? Ur 2 e 2 ^ \ 

e 2/3 

Ur 2 e 2 "< -r 2 e 2 "i 

\ -r 2 e- 2 t sin 2 9 J 

with U, V, (3, 7 being functions of u, r and 9 was first given in @. We shall be using the forms found 
for the GR metric functions, and all derived quantities, as given in Q. 

In order to find the NGT generalization of the metric tensor ( |l3| ) , we require that the symmetric 
part of the NGT metric tensor is formally the same as the GR metric tensor. We then impose 
the spacetime symmetries of the symmetric metric on the antisymmetric tensor g^y This is 
achieved by demanding that £g . 9[^ v ] = 0, where the Killing vector field £^ is obtained from 
£^ ^9(fiu) = 0. The solution to this equation for the metric (|l3|) yields the single Killing vector 
field, = = sin 2 0(90. Imposing £? g^] = yields: 

£(i)fyg[H + %3]&£(i) + 9[3v]d^(l) = • ( 14 ) 

This equation gives d^g^ = 0, but does not exclude any components of g\^ u }- This is markedly 
different from the static spherically symmetric case, where the above procedure excludes four of 
the six components of g\av\ ■ 

To make the problem tractable, we need to determine which g[^ u ] should be excluded. The 
criteria to be applied are that our solution must reduce to the usual NGT static spherically sym- 
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metric solution when passing through the equilibrium position, and that for each set to zero, a 
corresponding field equation must also vanish. To accomplish this, we note that the imposition of 
axi-symmetry splits the antisymmetric field equations ( 10b| ), Qllb| ) into two sets of three indepen- 



dent equations. (This can be seen directly from the block-diagonal form of the GR metric). The 
first set explicitly involves the three functions g^ ,g^ ,g il2 ^'- 

(V=gg [H ) =o (0 = 0,1,2), 



#[oi,2] = o. (15) 
These four equations are not independent owing to the existence of the one identity: 

^gg [H ) =o 0^ = 0,1,2). (16) 



The second set of four equations explicitly involves the three functions g^ , g^ , g^ : 







These four equations are also not independent due to the one identity: 



. 



(17) 



(18) 



L [3(U,H>P 

We note that eliminating one complete set of three functions simultaneously eliminates the three 
corresponding equations. However, eliminating one or two functions from either set fails to reduce 
the number of equations and leads to an over constrained system. Moreover, a Killing vector 
analysis performed at the instant when the system passes through its equilibrium position, shows 
that the set {g^ 01 * ,g^ ,g^} must reduce to g' 01 ', while the set {g^ ,g^ ,g<- 32 *} must reduce to 
gl 32 ]. Since we know from the exact spherically symmetric solution that asymptotic flatness requires 
g[ 32 l = 0, we must eliminate the entire second set of three functions {g^ , g [31] ,g [32] }- 
In view of the above, the NGT generalization of the metric tensor (fl3|) is: 

( y r ~l e 2f3 _ JJ 2 r 2 e 2 1 e 2? +UJ Ur 2 e 2 7 + X Q 

e 2f} -oj a 

Ur 2 e 2 ^-X -a -r 2 e 2 ^ 





9iw 



r 2 e 27 sin 2 1 



(19) 
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where to, A and a are functions of u, r and 9. The contravariant metric tensor is given to order h 2 
(where h = oj, A or a) by: 



where 



flU 





e -2/3 +g [W] e - 



-2(3 + g [m\ 
~W Vr -i e - 



,[20] 



J01] 



J02] 



,[12] 



-(o> — aU)e 



-4/3 



-2/3-27 -2 

-cxe M 'r , 



^[02] 

2 ^[7 + ot 12 ! 



_ e -2 7r -2 







_ r -2 g 27 g j n - 



- {Ur 2 {aU - io)e-^ + e- 2 ^{\ - aVr- 1 )) v 
det (g^u) 



-r 4 sin 2 0e 4/3 . 
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4 Solving the Field Equations 

We shall solve the field equations (|l^) using the GR background given by 

1 



2 -2 

--c r + ... 
4 



-(eg + 2ccot 0)r~ 2 + (27V + 3cc + 4c 2 cot 0)r -3 + 
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Ng + iV cot 9 - c/ - 4cc e cot 9 - ^c 2 (1 + 8 cot 2 i 



y=r-2M 

where c(u,9),N(u,9),M(u,9) are functions of integration and C(u,9) satisfies: 



r 1 + .. 



(21a) 

(21b) 
(21c) 
(21d) 



4C )W = 2c 2 c n + 2cM + N cot9 — Ng. 

We begin by solving the three field equations (V — 99^^ J = 0> if 1 = 0>1>2) with the skew 
functions expanded in inverse powers of r. To this end, it is convenient to work with the following 
linear combinations of oj, A and a: 



(uj - aU)e~ 213 , 
(A - aVr- l )e- 2 ^ 

27 
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(22a) 
(22b) 
(22c) 
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The field equations then become: 

(r/r 2 ) ir sin 9 + (5 sin B) fi = , (23a) 
5,u + K,r ~ (Ur 2 r]) tr = , (23b) 
r 2 r/ >n sin 9 + r 2 (Ur]sm 9)^ — (k sin 9) t g = . (23c) 

Since g^ v must transform under the Poincare group when r — ► oo, we require that ry, k and 5 have 
the following expansions in 1/r. 

t? = m/r + n 2 /r 2 + ... , (24a) 
<5 = d + d!/r + d 2 /r 2 + ... , (24b) 
« = fc + fci/^ + ^A 2 + ... • (24c) 

We begin by solving for the coefficients of /i at lowest order in 1/r. It is only when we find these 
coefficients that we can be sure that the next terms in the expansion of h are smaller than the h 2 
terms that we have neglected. 

Equation ( |23b| ) yields do, u = 0, and since do = when the system passes through its equil- 
librium position, do = always. Equation fl23a|) then gives n\ = 0. Continuing to next order, we 
find from equation ( |23a| ) that (d\ sin0) $ = 0. This equation for d\ demands d\ = in order to get 
a solution that is regular on the polar axis. The next equation in the set, ( p3b[ ), yields k\ = c?2,n 
while the equation ( |23c[ ) yields ri2, u = sin -1 9{ko sin^)^- 

To gain some understanding of the expansions thus far, it is instructive to calculate the NGT 
charge, L 2 , of the system 

L 2 = ^-f (V=gg l0u] ) d 3 x = - f* n 2 sin 9 d9 =< n 2 > , (25) 



4tt J V v JJ J, v 2 
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where the brackets <> denote the angular average. To maintain a notation consistent with other 
work in NGT we define n 2 = l 2 Q(u, 9) where I 2 is a constant (with respect to r, 0, u and 9) with the 
dimensions of a [length] 2 and is identified as the usual NGT charge, and < Q >= 1 for times when 
the system passes through its equilibrium position. Notice that L 2 U = by dint of the condition 
n 2 ,u = sin -1 9(ko sin 9) ; g. This expresses the important fact that the NGT charge of a body cannot 
be radiated away, and is analogous to the situation in electromagnetism where the electric charge 
of an antenna does not change. 
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We can now use the results of the lowest order expansion to find out how far we can safely 
expand h in inverse powers of r before needing to incorporate h 2 terms. This can easily be ac- 
complished by studying the exact form of g^ u and g. We find that the linear order can accurately 
determine 77, k and 5 as follows: 

rj = Ql 2 r- 2 + Rl 3 r- 3 + Sl 4 r~ 4 , (26) 

K = II + Jfr- 1 + Kl 3 r- 2 , (27) 

S = Al 3 r~ 2 + Bl 4 r~ 3 , (28) 

(Only even powers of I will appear in physical quantities such as the curvature tensor). The 
functions of integration are dimensionless functions of both u and 6. They satisfy: 

_ (£$L|k, (29a) 

suit) 

J = IA }U , (29b) 
„ (^sin^) a 

R = { . )fi , 29c 
smU 

„ (Jsm9) a 

K = !%L_(°* (29e) 

S = (29f) 

1 /sine (c, fl + 2ccotfl + 1^0 \ 

#,u = -t-t ^ • 29g) 

Notice that ( |29cj ) , ( |29dD can be combined to give ( |29b| ) and that ( |29| ) , (|29g| ) can be combined to 
give ( |29ej ) . This is a consequence of the identity ( |l6|) . 

Substituting the first terms in the above expansions into ( |12b| ) gives: 

3/ 

W m = -^{2lJ + 2cI + lQ, e ),u + 0(r- 3 ) , (30a) 

Wrcil = 0(r- 4 ), (30b) 



31 



W M = ^(2U + 2cI + lQ, e ) + 0(r- 4 ) . (30c) 
The identity W^^^y = then gives 

(2l 2 J + 2clI + l 2 Q fi ), u = 

2l 2 J + 2clI + l 2 Q, e = l 2 f{0) . (31) 
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The field strength F, iup is 



2 



F012 = + 0(r- 3 ) , (32) 



and since -F012 = when the system passes through its equilibrium position, f{6) = 0. 
Our results can be summarised as folows: 





= 0(r 


" 2 ) + 


WW] 


= 0(r 


- 3 ) + 




= 0(r 


- 3 ) + 



(33a) 
(33b) 
(33c) 



5 Calculating the Radiation Flux 

The rate of energy loss for an isolated body in NGT is given by 

-R 2 I t 0l hidn , (34) 



dE 
~dt 



where the integration is over a sphere of radius R in the wave zone, hi is an outward pointing unit 
vector, and t^ u is the energy-momentum pseudo-tensor. The skew contribution to t^ u is given by 

GD: 

A f imx0 f v __GR l iv F 2\ 
6 skew _ ^2 a/3 \<i y > 



+ -. 



2 



(2hW GR g u(3 W^ a] + 2h^ GR g^W [M - GR g^h^W M ) . (35) 



Inserting our expressions (|33c| ) , ( |33aD and (|33b| ) for F^a, and W M , we find that 



'!>,!/], we nna xnar r skew = 
0(r~ 5 ) so there is no skew contribution to the radiation flux in NGT. The only non-vanishing 
gravitational radiation flux in NGT comes from the leading order terms of the symmetric sector, 
which reproduce the usual GR radiation formula. 

As explained in the introduction, the linearized wave equations allow us to obtain the radiation 
pattern of any generic source by linearly superimposing a suitable combination of differently aligned 
axi-symmetric solutions. Since we know that the radiation flux is identical in NGT and GR for the 
axi-symmetric case, we know that the flux will be identical in general. This line of reasoning can 
be rigorously supported by generalizing the work of Sachs ||, who found that the axi-symmetric 
wave solution of Bondi, van der Burg and Metzner j|| contained all the essential features needed 
to describe gravitational radiation from a generic source of no particular symmetry. 
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6 Conclusions 



We have proved that the flux of gravitational radiation from an isolated source in NGT is identical 
to that found in GR, a result that clearly refutes the claims made by Damour, Deser and McCarthy. 

In many respects our result is not unexpected. The two conserved sources for an isolated 
system in NGT are L 2 , the NGT charge, and T^ u , the energy momentum tensor. Since I? has the 
dimensions of a [length] 2 , while the constant M in the GR sector has the dimensions of a length, 
we immediately expect the skew sector to be shorter ranged than the symmetric sector. Moreover, 
from the viewpoint of the phenomenological model where I? is taken to be proportional to the 
fermion number, we expect that L will be a constant since the fermion number of an isolated body 
is a constant. This was borne out in our solution. Thus, there will be no dimensionless quantity L 2 U 
available in the skew sector unlike the symmetric sector where the dimensionless quantity M u is 
available. Using dimensional analysis, we can then predict the form of the Riemann tensor, which 
must have the dimensions of [length] -2 and contains at most two time derivatives. The symmetric 
sector can contribute super-radiative 1/r terms and radiative 1/r 2 terms, while the skew sector can 
only contribute non-radiative 1/r 4 (and lower) terms. 
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